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The  apparent  regularity  in  energy  distribution  among  Internal-wave  apace  end  time  scalea 
may  be  a saturation  phenomer.in.  This  paper  develops  a model  of  the  saturated  internal-wave  field 
baaed  on  the  hypothesis  thst  the  saturated  quantity  is  fiz/N  , the  locally  normalized  squared  vor- 
ticity.  This  quantity,  which  la  a generalized  inverse  Richardsot  number  for  spatially  varying 
flow,  is  shown  in  the  linearized  approximation  to  be  the  dc.  sity  cf  a conserved  "excitation"  1; 
the  excitation  la  asstxned  to  be  uniformly  partitioned  among  the  novmsl  modes  in  a saturated  field. 

The  resulting  energy  distribution  differs  from  that  of  Garrett  and  Hunk  in  that  energy  la 
assigned  preferentially  to  the  lowest  modes.  Among  the  properties  derived  for  this  distribution, 
and  subject  to  verification  by  ocean  measurements,  is  a universal  law,  independent  of  the  strati- 
fication profile  N(x) , for  coherence  spectra  of  iaopycnala  in  the  jpatial  domain. 
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1 . INTRODUCTION 

'V  /wvwvv/vAi 

Evidence  is  accumulating  that  the  distribution  of  ocean  internal 
wave  energy  among  space-time  scales  may  follow  a universal  law.  Garrett 
and  Munk  (1972.  1975)  have  proposed  tentative  forms  for  this  distribution 
that  agree  fairly  well  with  a variety  of  spatial  and  temporal  measure- 
ments taken  in  several  locations.  It  is  tempting  to  ascribe  such  a uni- 
versal distribution,  if  it  indeed  exists,  to  a saturation  phenomenon,  in 
w^ich  energy  in  excess  of  the  saturated  value  at  a given  wavelength  and 
frequency  is  quickly  dissipated  through  nonlinear  exchange  or  turbulent 
breaking. 

In  another  paper,  Garrett  and  Munk  (1972  A)  explore  this  idea 
further.  They  argue  that  a saturated  distribution  will  be  marginally 
stable,  experiencing  instabilities  only  occasionally  to  balance  the  slow 
input  of  energy.  Borrowing  the  inverse  Richardson  number,  Ri  \ from 
the  theory  of  parallel  stratified  flow  as  a measure  of  instability,  they 
find  that  the  values  implied  by  their  distribution  are  of  an  order  that 
would  indeed  suggest  marginal  stability. 

Suppose  we  take  this  idea  a step  further:  If  the  limiting  distri- 

bution is  determined  by  physical  saturation  such  that  Ri  , or  some  other 
suitable  measure  of  excitation,  attains  a uniform,  universal  statistical 
level,  we  ought  to  be  able  to  use  this  information  to  deduce  a priori 
the.  form  of  the  distribution.  In  this  paper,  such  a saturation  hypo- 
thesis is  proposed  and  its  consequences  worked  out  in  detail. 

For  reasons  that  will  be  explained,  a modified  measure  of  excita- 
tion better  suited  to  internal  waves,  the  vortex  intensity , is  introduced 
in  place  of  the  traditional  Ri  The  vortex  intensity  resembles  Ri  ^ 
except  that  total  vor "lefty  is  substituted  for  velocity  shear.  As  will 
be  seen,  the  assumption  of  a statistically  uniform  level  of  excitation 
throughout  the  thermocline  leads  to  an  energy  distribution  that  is 
physically  reasonable  and  not  in  obvious  contradiction  to  existing  spatial 
and  temporal  measurements,  yet  quite  different  in  certain  respects  from 
the  distribution  of  Garrett  and  Munk.  The  major  difi«.rence  is  in  the 
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partitioning  of  energy  among  the  normal  modes  or  vertical  wavenumbers, 
and  this  difference  is  not  easily  distinguishable  observationally  except 
through  cross-spectral  analysis  of  multi-depth  data.  Curiously,  the 
model  of  uniform  vortex  intensity  developed  here  yields  exact  analytic 
predictions  of  spatial  autospectra  and  vertical  coherence,  independent 
of  the  density-stratification  profile,  and  on  this  account  the  model  is 
very  directly  testable. 

The  treatment  given  here  assumes  no  steady  shear  current  and  no 
rotaticn  (no  inertial  waves).  One  horizontal  dimension,  x,  is  used  for 
notational  convenience,  but  the  generalization  to  two  dimensions  is 
straightforward . 


» 


( 


i 


I 


2. 


LINEARIZED  NORMAL  MODE  EQUATIONS 

The  field  c(x,z,t)  of  vertical  fluid  displacement  is  decomposed 
in  the  usual  way  into  normal-mode  oscillations  of  horizontal  wavenumber 
k,  phase  speed  c,  and  frequency  m - ck, 

;(x,z,t)  - am(V.)  $ra(z)exp[ik(x  T c^t)],  (1) 

k,m 

the  eigenfunctions  ^ satisfying 


+ (c“2  N2(z)  - k2)  <j>  = 0 (2) 

,2m  is 

dz 

and  vanishing  at  the  top  and  bottom  of  the  water  column  (Phillips  1966, 

Ch.  5).  It  is  convenient  to  regard  k as  an  arbitrary  but  fixed  separa- 

_2 

tion  parameter  on  which  the  discrete  eigenvalues  cm  and  eigenfunctions 

ij>  depend;  the  eigenfunctions  then  make  up  an  orthogonal  and  complete 
m 

set.  We  use  the  normalizing  convention 

L n24>  dz  * 6 , (3) 

m n mn 

according  to  which  the  completeness  relation  is  symbolically 

oo 

N2(z)  V''  <J>  (z)«f>  (z’)  = 6(z-z’).  (A) 

m m 

m=l 


1 


i 


Tht  associated  fields  of  horizontal  velocity  u and  vertical  velocity  v 
are,  for  each  mode, 
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(5) 


r J 


u ■ ± a c <f>' 
u m m m 


v * T i a k c $ , 
m m m m 


where  <f>'  is  shorthand  for  dd>/dz  and  where  the  factor  exp[ik(x*c  t)]  has 

m J 

been  absorbed  into  the  mode  amplitude  a for  brevity. 

We  assu„i»*  -hat  the  internal  wave  fields  found  in  nature  are  homo- 


geneous and  random  so  that  their  statistical  properties  are  adequately 

2 

described  by  the  me.'  .-square  mode  amplitudes,  <laml  >»  or  more  precisely 


by  spectral  densities  A defined  over  all  wavenumbers, 

m * 


<|a  > A (k)dk. 

id 


Our  aim  is  to  determine  these  amplitude  spectra,  on  the  one  hand  from 
a principle  of  saturation,  and  on  the  other  hand  from  empirical  meas- 
urements of  assoc  Uted  displacement,  current,  and  energy  spectra  to  which 
the  amplitude  spectra  are  uniquely  related. 

For  example,  the  total  wave  energy  in  a water  column  of  unit  sur- 
face area,  at  a given  k. 


y fin  + v2  + N2£2]dz, 


(6) 


becomes  with  the  aid  of  Equations  (1)  through  (5) 

£ W* 


E = p 


(7) 


m 


so  that  the  spectral  density  of  energy  in  each  mode  is 


E (k)  - p A (k), 

in 


(8) 


VORTEX  AMPLITUDE  AND  INTENSITY 


To  treat  the  topic  of  internal  wave  excitation  as  generally  as 
possible,  we  need  an  absolute  dimensionless  measure  of  local  amplitude 
that  plays  a role  for  internal  waves  analogous  to  that  of  slope  for  sur- 
face gravity  waves.  One  candidate  is  the  ratio  of  current  shear  to  local 
Vai sal a- Brunt  frequency,  u’/N(z),  whose  square,  the  inverse  Richardson 
number  Ri"  , is  the  measure  of  intensity  used  by  Garrett  end  Munk  to 
assess  the  degree  of  saturation  of  their  proposed  energy  distribution 
(1972  A)  and  earlier  by  Phillips  (1966,  Ch.  5)  to  estimate  an  upper 
unit  to  the  possible  excitation  of  the  first  (interface)  mode.  This 
normalized  amplitude  has  a serious  defect  when  applied  to  waves  in  a 
vertical  continuum,  as  can  be  seen  in  its  normal  mode  expansion  via 
Eqs.  (2)  and  (5),  shown  here  at  a single  wavenumber  k for  simplicity: 

u' (x,z,t)/N(z)  - N_1 (z) 


- N_1(z) 


This  amplitude  becomes  infinite  in  regions  of  vanishing  stratification 
(except  at  k - 0),  regions  where  the  local  flow  is  stable  potential  flow. 
The  infinities  arise  because  of  the  k2  terms  in  the  above  expression,  a 
clue  that  the  shear  amplitude,  borrowed  from  stability  theory  for  parallel 

shear  flow,  has  not  been  properly  adapted  to  flows  having  horizontal 
periodicity. 

This  difficulty  vanishes  when  the  shear  u’  is  replaced  by  total 
wave  vorticity 


a c «p 

m mm 


m 


(9> 


a 


9v 

9x 


9u 

9z 
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to  define  a vortex  amplitude 


e E fl/N; 


(10) 


the  expansion  of  this  quantity  is 


e(x,z,t) 


N_1(z) 


mm  m id 


" N(Z)  J^Vm^m’  (ID 

m 

showing  that  the  amplitude  is  well  defined  regardless  of  the  detailed 
shape  of  the  stratification  profile.  The  generalized  definition  of 
local  intensity  is  then 


e2(x,z,t)  - n2/S2,  (12) 

which  reduces  to  Ri  in  the  long-wave  limit,  k^O. 

Equation  (11)  above  indicates  that  the  vortex  amplitude  and  inten- 
sity are  confined  to  stratified  layers  even  though  part  of  the  kinetic 
energy  may  reside  elsewhere  in  unstratified  regions.  Intuitively,  the 
question  of  instability  comes  down  to  the  mutual  interaction  of  vortex 
layers,  an  interaction  in  which  neighboring  unstratified  flow  can  be 
viewed  as  playing  a more  or  less  passive  role.  On  this  basis,  the  vor- 
tex intensity  would  seem  to  be  a plausible  measure  of  the  local  tendency 
to  instability. 

As  a general  descriptor  of  internal  wave  excitation,  the  vortex 
intensity  has  another  very  desirable  property; 
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Vortex  intensity  is  a conserved  dencity. 


Consider  a system  of  internal  waves,  again  at  fixed  k for  simpli- 
city, with  an  arbitrary,  time-dependent  vertical  structure.  We  can  define 

2 

a total  excitation  I per  unit  surface  area  by  integrating  e over  the 
vertical  column, 


I 


(13) 


a definition  analogous  to  that  for  total  energy,  with  e replacing  dis- 
placement and  velocity.  Squaring  Eq.  (11),  we  get 


I 


/£ 

m,n 


* -1-1  2 
a a_c  c_  $ N $ dz. 


m n m 


m 


(14) 


and  by  the  orthogonality  rule  of  Eq.  (3)  we  find 


I 


-2 

c * const: 

m 


(15) 


In  the  limit  of  linearized  dynamics,  the  excitation  I belonging  to  a 
group  of  waves  remains  unchanged  as  the  group  propagates  vertically 
through  different,  possibly  dissimilar,  parts  of  the  thermocllne. 

The  excitation  has  dimensions  of  length  denoting  both  the 
level  of  intensity  and  the  total  depth  of  the  water  column  involved. 
Equation  (15)  suggests  that  the  excitation  is  partitioned  among  the  modes 
according  to 


I - 

m 


(16) 


which  together  with  Eq.  (7)  implies  that  the  mode  energy  and  excitation 
are  related  by 
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! 


2 _ 

pc  I . 
m ra 


PROPOSED  ENERGY  PARTITIONING 

We  oow  pass  to  a statistical  description  and  replace  a given  sys- 
tem of  waves  by  an  ensemble  of  systems  in  which  the  mode  amplitudes  are 
random  in  phase  and  amplitude.  The  ensemble  will  be  stationary  in  time 
(at  least  with  respect  to  mean  squares  and  products  such  as  energy, 
excitation,  and  space-time  covariances)  if  the  mode  amplitudes  are  uncor- 
related: 


* i 1 2 

<a  a > * < a > 6 
m n m ' mn 


Defined  thus  in  terms  of  mean  square  amplitudes,  the  ensemble  is  formally 
in  equilibrium  under  linearized  dynamics.  This  should  be  distinguished 
from  physical  equilibrium,  which  is  attained  when  the  mean  square  ampli- 
tudes slowly  adjust  themselves  to  bring  the  generation  mechanisms,  non- 
linear energy  exchange  among  modes  and  wavenumbers,  and  dissipation  into 
balance.  We  assume  that  the  time  scale  for  the  processes  of  physical 

equilibrium  is  much  longer  than  the  oscillation  periods  tu  so  that  the 

m 

stationary  ensemble  with  linearized  dynamics  is  a useful  model  of 
equilibrium  over  many  periods. 

Can  we  bypass  the  intervening  physics  and  directly  characterize 
the  saturated,  "fully  aroused"  internal  wave  field  by  means  of  the  vor- 
i ex  amplitude?  We  have  argued  that  the  vortex  intensity  is  a plausible 
measure  of  local  excitation  relative  to  instability,  and  have  shown  that 
it  is  the  density  of  an  invariant,  the  so-called  excitation  1.  A rea- 
sonable guess  is  that  in  a saturated  field,  the  excitation  is  uniformly 
distributed  in  depth  and  among  length  scales.  We  make  the  following 
hypothesis: 
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At  a given  wavenumber,  the  excitation  is  partitioned  uniformly  among  the 
modes . 


This  proposed  distribution  implies  the  following  relation  among 
modal  mean-square  amplitudes  and  energies: 


<1  > - I , » - 1,2, 

m o 


(19a) 


<|a  |2>  - c2  I , 
' m'  mo 


(19b) 


<E  > * pc  I , 

m mo 


(19c) 


where  I is  the  universal  excitation  per  mode. 


Totcl  energy  is  finite. 


The  total  mean  energy  is 


<E>  = pi 


2 v 


now  c^  is  a declining  function  of  m (see  Figure  1),  which  in  the  limit 
of  high  mode  numbers  is  not  badly  approximated  by  the  WKB  values 


c ca  n 1(m  - ^ 1 / N(z)dz. 

m L J 


2 -2 

Thus  c^  ~ m , and  even  for  infinitely  many  excited  modes,  the  total 
energy  is  finite. 

Practically  speaking,  some  sort  of  mode  cutoff  is  required  to  keep 
the  total  excitation  from  diverging.  Viscosity,  if  nothing  else,  will 
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2 

damp  modes  of  vertical  scale  h where  id  h is  lass  than  the  kinematic 

m 


viscosity.  It  is  important  to  observe  nevertheless  that  the  total 
energy,  and  other  properties  such  as  predicted  covariances  of  displace- 
ment and  velocity,  will  be  very  weakly  dependent  on  mode  cutoff,  so  that 
for  most  purposes  the  cutoff  can  be  ignored. 

The  vox 'ex  amplitude  is  statistically  unifoim  in  depth  and  vertical  wave- 
number. 

The  two-depth  covariance  of  the  random  vortex  a lplitude  field 
e(x,z,t)  has  a simple  form,  as  deduced  from  expression  (11)  via  the  com- 
pleteness relation: 


I 


<e(x,z,t)e(x,z' ,t)> 


N 


E*  -1  -1 

<a  a >c  c <j>  (z)$  (z') 
mn  m n m rn 

m,n 


^r*ffiu)*m(z’) 


- Iq  6 (z-z ' ) . (22) 

Amplitudes  at  different  depths  are  completely  uncorrelated,  and  the  field 
e is  seen  to  be  spectrally  white  with  respect  to  vertical  wavenumber: 

The  vortex  amplitudes  are  in  a state  of  "maximum  chaos". 

DISTRIBUTION  OVER  HORIZONTAL  SCALES 

To  generalize  the  foregoing  from  a particular  horizontal  wave- 
number  to  all  wavenumbers  simultaneously,  one  merely  replaces  every  mean 
square  quantity  by  its  corresponding  spectral  density.  The  f ndamental 
spectrum  is  then,  by  hypothesis,  the  universal  excitation  spectrum  I(k) 
and  the  related  mode-amplitude  and  energy  spectra  are 
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3.  Some  Consequences  of  the  Proposed  Distribution 

VWVWVlAWiA/lAAA/WWWi/W\MwiWl/WWWVi/W\j 


A uniformly  excited  internal-wave  field  will  have  certain  observ- 
able characteristics  that  make  it  easy  to  distinguish  by  means  of  the 
appropriate  measurements.  As  one  might  guess,  the  most  definitive  meas- 
surements  will  be  these  capable  of  directly  resolving  horizontal  wave- 
number,  i.e.,  multi-depth  towed  measurements,  since  the  proposed  distri- 
bution assigns  energy  over  spatial  scales.  Temporal  data  from  moored  or 
floating  multi-depth  measurements  can  provide  additional  checks,  but  as 
will  be  seen,  probably  cannot  by  themselves  confirm  the  distributions. 

A striking  and  unexpected  property  of  the  uniformly  excited  field 
is  the  existence  of  simple  closed-form  predictions,  independent  of  N(z), 
for  spatial  autospectra  and  vertical  coherence  spectra  for  the  quantity 
C,  vertical  displacement. 


SPATIAL  AUTO-  AND  CROSS-SPECTRA 


The  Instantaneous  values  of  a random  displacement  field  <;  can  be 
expressed  formally  as  an  integral  over  wavenumber  of  the  stochastic 


amplitudes  a (k) , 

ID 


C(x,z) 


»<kH» 


(z,k)exp[ikx]dk, 


where  the  amplitudes  have  spectral  densities  A (k)  and  are  mutually  uncor- 


related at  different  k and  m: 


<a  (k)a  (k')>  * A (k)6  6(k-k') 

m n ~m  mn 


It  follows  that  the  cross  spectrum  of  Q at  two  depths,  z ^ ana  z 2,  is 


given  by 
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(27) 


~i2,k) ' £4,(k)W*„<22> 

m 


with  the  single-depth  autospectrum  represented  by  the 

2 2 ’ zi; 


special 


case 


zi(k)  ' £4»<k,*m<2i: 


(28) 


m 


By  sq.  .23)  for  umiform  erclt.tlon.  ^(k)  - c2(k)l(k),  we  then  have 


^(k)  ■ i(k>  D <ww- 


(29) 


m 


This  is  an  equation  for  the  (symmetric)  dependence  of  . „„  the  t„„ 

arguments  z and  z„  at  each  k Tf  „ k ~12 

twice  „,tl  2 8 h°ld  Z2  “1  differentiate 

twice  with  respect  to  2j,  « g.t  vU  Eq.  (2/ 


_3 

dz 


i. 

hi- 1 Y 

m 


and  subtracting  k2  times  Z 

~12’ 


(~ 


9 2 ’ k2  ) l 


~12  " i " (ll>  X)m(z2>*m<zl> 


m 


' I «(z1-z2). 


(30) 
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Now  since  £ vanishes  at  the  top  and  bottom  of  the  water  column,  the  same 
must  hold  for  80  that  Eq.  (30)  has  the  unique  solution 

m lsinh  kz<exp(-ks>)  (31) 

where  z<t  z>  are  the  lesser  and  greater  of  the  depths  z ± and  z2.+  An 
equivalent  expression  in  terms  of  depth  separation  bz  * z -z  is 

~12 (k)  " 1 I(k)k  1 exP (-kAz) [l-exp(-2kz<) ] . (32) 


Autospectra 


The  displacement  autospectra  are  a one-parameter  family 

^l(z)  " 2 £<k)k  1[l~exp(-2kz)]  (33) 

directly  related  to  I(k).  The  asymptotic  limits 

| zl(k),  k < z"1 


2 k £(k) » k > z , (34) 

indicate  an  increase  in  negative  spectral  slope  of  unity  from  low  to 
high  wavenumber,  with  the  transition  occurring  around  k ~ z~ The 
spectral  densities  collapse  to  a single  value  for  depths  greater  than 
k , and  are  proportional  to  z for  depths  smaller  than  k"1.  This  behavior 


^The  form  shown  is  for  an  effectively  infinite  bottom  depth  d,  i.e.,  when 
kd>>l.  Otherwise,  replace  exp(-kZ>)  by  sinh  k(d-Z>)/sinh  kd. 
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is  illustrated  in  Figure  2,  for  an  assumed  excitation  spectrum  of  the 
form  I(k)  ~ k 1,S.  The  tow  spectia  of  Charnock  (1965),  well  known 
through  their  reproduction  *n  the  monograph  of  Phillips  (1966),  are 
accounted  for  nicely  by  the  prediction  above  (uee  Figure  3)  if  we  allow 
ourselves  the  freedom  to  adjust  the  level  and  spectral  slope  (-1.5)  of 

I(k).  Of  course,  this  i.5  barelv  more  than  a consistency  check,  but  we 
do  observe  a slope  increase  at  kz  ~ 1.0. 

Vertical  coherence. 


The  two-depth  coherer ce  spectra  predicted  by  Eqs.  (31)  and  (33) 

are 


~12V 


(k)  i 


(£i 


1/2 


__l/2 


exp(-kAz) 


sinh  kz< 
sinh  kz> 


(35) 


In  the  low  and  high  wavenumber  limit,  the  behavior  is  approximately 

( <VZ>>1/2.  k < z^1 

£i2(k)- 

( exp(-kAz),  k > z<  \ (36) 


indicating  moderately  high  coherences  at  low  wavenumber  when  the  depths 
are  not  too  different  and  a rapid  decline  in  coherence  around  kAz  ~ 1. 
The  horizontal  wavelength  for  p • 0.5  is  about  9(Az).  See  Figure  4. 


Transfer  spectrum. 


A quantity  usually  used  in  other  applications,  the  transfer 
spectrum 
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DISPLACEMENT  SPECTRAL  DENSITY.  M2HR/CYCLE 


Figure  3 


TOW  DIRECTION, 
DEGREES 

• 068,  248 
‘113,  293 

* 158,  338 
° 205,  023 


the  reported  depth  and  tow  speed  of  75  m and  11  5 kt  iiZL  a at 
"Ution  U Independent  ofV),  the  rCntUe  Ldli 

75  m^enth^f^l6  °f  f]9ljre  1 ’ whose  maximum  Is  close  to  the 
75  I*CSrn0ckJf  measurements,  the  first  mode  accounts 

for  nost  of  the  predicted  displacements. 
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KZ,,  RADIANS 


Figure  4.  Spectra  tf  Coherence, for  Spatial  Series  of  Vertical  Displacement 
at  Two  Depths  In  a Uniformly  Excited  Field.  These  spectra  are  a 
universal  function  of  KZ,  ai  d Z2/Z1 . 


L9 


Il2(k)  1 212(ky41<k)»  (37) 

has  a particularly  simple  predicted  f*r-m  ( z 1 is  the  lesser  depth): 

T12(k)  - exp(-kAz)  . (38) 

The  behavior  of  the  multi-depth  autospectra  and  transfer  spectra  pre- 
dicted above  will  make  a uniformly  excited  internal-wave  field  easy  to 
recognize  on  the  basis  of  fairly  simple  analysis  of  multi-depth  towed 
measurements,  [a  towed  thermistor-chain  experiment  has  been  conducted 
very  recently  (Nelson,  1975)  and,  for  the  first  time,  data  of  sufficient 

scope  and  quality  for  analysis  in  both  horizontal  and  vertical  dimensions 
will  shortly  be  available.] 

Treatment  of  towed  measurements 

The  spectral  densities  treated  above,  which  are  intended  to  repre- 
sent distributions  over  scalar  wavenumber  in  two  dimensions,  are  not  the 
same  as  tow  spectra  obtained  by  Fourier  transform  and  should  not  be 
strictly  compared.  The  expected  covariance  of  isotherm  levels  at  two 
depths  zr  z2  and  two  horizontal  positions  x^  x2  along  the  tow. 


C12(VX2>  H <^^x1»z1)?(x2»z2)>»  (39) 

is  related  to  the  (equivalent  one-dimensional)  cross-spectral  density 
Z1200  by 

C12  (xi"x2)  " 27 //*i2<k>  exP[ik(x1-x2)cos6]dkd0  (40) 

where  kcos0  is  the  projected  periodicity  along  the  tow  direction  of  com 
ponents  making  an  angle  0 with  the  tow  direction-  Thus, 
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ci2(xrx2) 


(41) 


oo 

J Z12(k)  JQ(k|x1-x2 |)dk. 


and  the  proper  inverse  relation  is  the  Hankel  transform 


OO 

/ 


Zi2(k)  - k J C12(x)JQ(kx)xdx. 


(42) 


This  indicates  how  true  scalar  wavenumber  densities  can  be  recovered 
from  empirical  linear  covariances,  either  under  the  assumption  of 
isotropy  or  after  averaging  over  several  tow  directions. 


TEMPORAL  SPECTRA  AND  COHERENCES 


The  temporal  spectrum  S(oi)  of  a quantity  or  pair  of  quantities 
can  be  directly  inferred  from  the  spatial  spectra  Sm(k)  under  the  assump- 
tion that  the  modal  oscillation  frequencies  are  reasonably  well  repre- 
sented by  the  dispersion  relations  ^(k)  predicted  by  the  eigenvalue 
equation  (2).  Figure  5 is  a representative  dispersion  plot  for  the 
sample  N(z)  profile  of  Figure  1.  Since  the  variance  or  covariance  in 
the  interval  dio  is  a sum  of  contributions  from  all  modes  in  the  propor- 
tion 


m 


where  ^(k^)  ■ uj,  we  have 


S(o>) 


ES  (k  ) 

~m  m 


dk 


-1 

k 

m 


(43) 
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K,  CYCLES/KM 


Figure  5.  Modal  Dispersion  Frequencies,  «m(k)»  For  the  Profile  of  Figure  1.  These  were 
obtained,  along  with  the  associated  eigenfunctions  ) , by  numerical 
solution  of  the  eigenvalue  equation  (2),  and  used  in  the  sample 
calculations  of  temporal  spectra. 


For  displacement  (i.e.,  isopycnal)  cross-spectra  at  iwo  depths,  the 
appropriate  expression  is 


~12(a))  * ^4(km)cg^  ^m(zl)^m(z2)  (44a) 

m 

Where  Cgm  is  the  modal  8rouP  sPeed  dc^/dk,  and  where  Am  and  <J>  are 
evaluated  at  k = ^(w).  In  a modal  sum  at  constant  frequency,  it  seems 
preferable  to  use  the  rescaled  eigenfunctions  <p  (z,u>)  = (c  /c  )^2  &,  (z  k 
which  comprise  an  orthonormal  set  for  the  time  domain,  as  shown  in  the 
appendix.  The  sum  then  reads 

^12(U))  =S^(km)Cmlv,m(zl)%1(z2) 

m 

’“EC  i(km)<Pm(zi)%1(z2)'  (44b) 

m 

The  autospectra  and  coherence  spectra  are  accordingly 

•M">  - “ E ^ <45> 


and 


£i2(uj) 


|z12(»)|/Cz1(»)z2(«)] 


1/2 


(46) 


Unlike  the  spatial  spectra,  these  quantities  depend  in  a subtle  way  on 
the  stratification  profile  implicitly  through  k^w),  and  no  general 
formulae  are  immediately  evident. 
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Certain  qualitative  features  can  be  inferred.  At  low  frequencies 


c — ► const, 
m 


TE 


const. 


so  that  if  the  excitation  spectrum  has  the  form  of  a power  law,  1 ~ k P, 
the  mode  coefficients  in  the  sums  (44,  45)  behave  as 


-p  p+1 

I c ~ w c ; 
~m  m m 


n 


.1 


thus,  in  the  low-frequency  limit,  the  auto-  and  cross-spectra  have  a 
power-law  dependence  with  the  same  slope  as  the  spatial  autospectra  at 
low  wavenumber.  Note  also  that  for  any  appreciable  spectral  slope,  say 
p ~ 1.5-2,  the  factor  cP+1  diminishes  rapidly  with  mode  order,  since 

ID 

c ~ (m-l/2)-i,  so  that  the  sums  would  tend  to  be  dominated  by  the  first 
m 

mode.  This  would  suggest  that  the  predicted  coherences  remain  appre- 
ciable. 

For  the  thermocline  of  Figure  1,  modal  eigenfunctions  have  been 
computed  numerically  and  used  to  construct  displacement  spectra  and 
coherence  spectra  corresponding  to  I(k)  ~ k ^ These  are  shown  in 
Figure  6.  The  autospectra  resemble  some  observations  in  that  they  show 
the  expected  slope  and  -he  cutoff  at  the  local  N(z)j  however,  because 
of  the  first-mode  dominance,  Lhe  low-frequency  values  scale  as  z rather 
than  as  N_1  at  shallow  depth,  contradicting  the  rule  derived  from  the 
WKB  approximation. 

The  predicted  dominance  of  the  first  mode  in  temporal  spectra 
means  that  temporal  data  are  likely  to  be  ineffective  in  determining 
contributions  from  modes  higher  than  the  first,  if  the  distribution  is 
truly  uniformly  excited.  Thus,  temporal  data  probably  cannot  confirm 
such  a distribution.  On  the  other  hand,  certain  serious  departures  from 
uniform  excitation,  such  as  a "top  hat"  distribution  of  the  kind  used  by 
Garrett  and  Munk;  c^n  be  detected  in  temporal  data. 
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NORMAL  MODES  OR  EQUIVALENT  CONTINUUM? 

Whether  one  discusses  energy  partitioning  in  terms  of  normal 
modes  or  an  equivalent  continuum  (m,  k)  is  partly  a matter  of  indivi- 
dual preference.  Some  argue  that  for  energy  at  large  vertical  wave- 
numbers,  the  discrete  normal-mode  picture  is,  at  best,  a nuisance,  and 
at  worst,  misleading.  To  the  extent  that  the  effects  of  shear  currents 
are  important,  those  normal-mode  components  whose  phase  speed  is  not 
substantially  greater  than  the  velocity  shear  across  the  thermocline  are 
inaccurately  represented.  However,  in  the  vanishing-shear  approximation, 
the  normal  modes  are  as  complete  a description  as  any,  and  they  have  the 
advantage  of  being  an  orthogonal  set  of  "degrees  of  freedom,"  independent 
oscillators  whose  excitation  can  be  systematically  compared.  In  retro- 
spect they  seem  to  be  essential  to  a rigorous  description  of  a "uniformly 
excited"  field  because  this  field  has  most  of  its  energy  in  the  lowest 
few  modes.  For  example,  at  fixed  k the  first  mode  has  as  much  as  75% 
of  the  total  energy,  and  at  fixed  u even  more. 

What  is  the  equivalent  continuum  energy  density  for  the  uniformly 
excited  field?  If  we  define  D(a),k)dmdk  as  the  energy  in  the  entire 

vertical  water  column  and  in  the  frequency  intervals  dui,dk,  then  D is 

related  to  E by 
m 


D(Vk)  E„<k)  (47) 

when  is  the  modal  frequency  spacing  at  fixed  k.  According 

to  the  WKB  approximation,  which  Is  Increasingly  trustworthy  at  higher 
mode  numbers,  the  eigenfunctions  are  determined  by 

k /<"2(Z>/V  1>1/2  dr  * (m  - i)„.  (48) 

N>u) 


Preceding  page  blank 
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dz  — n » 


To  estimate  Aw  , we  can  differentiate  by  m, 


,.,d  solving  for  Aw"1,  the  "density  of  states"  at  the  frequency  u,  we  get 

Aw"1  a (irc2^'1  /N(l-wV)'1/2  dz,  (49) 

m m I m 

2 

which,  combined  with  * pcmi’  yields 

D(w,k)  = pn  1k  1I(k)  j N(l-w2/W2)  ^ dz.  (50) 

The  continuum  density  obviously  depends  on  the  detailed  shape  of  the 
Vaisala- Brunt  profile,  N(z),  but  if  this  profile  is  not  pathological, 
the  integral  in  Eq.  (50)  depends  only  weakly  on  w,  so  that  at  fixed  k, 

D is  at  least  roughly  independent  of  w for  frequencies  distinctly  less 
than  the  maximum.  Further,  D approaches  a well-defined  limit  as  w-^0. 

COMPARISON  WITH  THE  DISTRIBUTION  OF  GARRETT  AND  MUNK 

For  simplicity,  we  limit  the  comparison  to  frequencies  w well 

below  the  maximum  and  set  the  inertial  i quency  to  zero  in  the  Garrett 

Munk  expression,  which  allows  us  to  assume  that  the  dispersion  curves 

w (k)  are  nearly  straight  lines: 
m 

r 


w (k)  - c k, 
m m 


1 


and  we  adjust  I(k)  ~ k to  predict  the  same  slope  value  for  isotherm 
spatial  and  temporal  spectra  as  used  by  Garrett  and  Munk.  We  can  then 
compare  continuum  density,  along  with  equivalent  modal  energy  an^  excita- 
tion. A third  distribution,  one  in  which  the  energy  is  partitioned 
equally  among  the  modes,  with  properties  intermediate  to  the  other  two, 
has  been  added  to  the  comparison  below  for  purposes  of  discussion: 


| i 


I 


Uniform 

Excitation 

D(m,k) 

k"3 

k2E  (k) 

k2i  (k) 
~m 

const. 

m 

<ao 

Energy 

Garrett- 

F.quipartition 

Munk  *72 

, -1  -2 

-3 

k ui 

0) 

const. 

(m-y) 

(m-^) 

VI 

-r-j 

VI 

Both  the  equal-energy  and  Garrett-Munk  distributions  require  a mode 
cutoff  j,  which  limits  excitation  in  the  (u,k)  plane  to  the  region 

< c*k  < (j-|)ui. 

Note  that  both  energy  and  excitation  increase  with  mode  number  in  the 
Garrett-Munk  distribution,  the  latter  quantity  very  rapidly. 

All  three  stributions  have  been  properly  adjusted  to  exhibit  the 
same  autospectra  for  both  stationary  and  towed  isotherm  data.  Fortu- 
nately, they  can  be  expected  to  predict  very  different  vertical  coherences. 
Foi  example,  compare  the  spatial  cross  spectra  predicted  for  the  uniform- 
excitation  model, 

~12^)  = * sinh  kz<  exp(-kz>). 
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with  those  predicted  by  energy  equipartition  Q^OO  “ £00) 

hi  " pl£(k)  ^ ♦»(*1Hb(z2)* 


ID“1 


These  latter  become  very  abrupt  functions  of  z^”z2  *n  high 

mode  cutoff  j, 


Zi2(k)-^p"1£(k)N(zi)"26(zi-z2). 


suggesting  an  absence  of  measurable  coherence  between  any  but  the  nearest 
isotherms.  This  comparison  is  a particular  example  of  how  predicted 
coherences  diminish  with  increasing  "modal  bandwidth",  and  indicates  how 
the  moderate  bandwidth  of  the  uniform  excitation  model  could  be  distin- 
guished from  the  higher  bandwidth  of  a Garrett-Munk.  distribution  in  multi- 
depth data. 
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EIGENFUNCTIONS  AT  CONSTANT  u 

Measurements  of  internal  waves  in  the  time  domain  are  more  con- 
veniently described  by  a set  of  eigenfunctions  belonging  to  a common 
frequency  ui.  When  ui  is  used  as  an  independent  parameter,  the  eigenvalue 
equation  reads 


II 


* 

m 


+ c'2[N2(z)  - 
in 


0, 


(A-l) 


which  implies  that  the  eigenfunctions  obey  the  orthogonality  condition 


/■ 


(N2-w2)<P 
m m 


dz  - 6 


mn 


(A-2) 


Because  of  the  normalization  chosen  at  m - n,  these  functions  are  related 
to  the  previously  defined  eigenfunctions  by  some  non-unit  constant  of 
proportionality. 


(z,u))  - y (w)*  (z,k  (w)),  (A- 3) 

m m m m 


where  k (w)  lies  on  the  dispersion  curve 
m 


k (u>)  = w/c  (u>) . 
m m 


(A- 4) 


(See  Figure  5.)  At  constant  frequency  the  eigenfunctions  are  complete 

2 2 

only  on  the  portions  of:  the  vertical  domain  where  N -uj  >0;  formally. 


[n2(z)-cj2]  <p  (z)x(z')  = 6(z-z').  (A-5) 

m m 

m*l 


, 


| 


A-l 


The  proportionality  constant  Ym  niay  be  related  to  the  group  speeds. 


c = du)  /dk 

gm  m 


c [l  + -r  dc  /dk] 
m cm 
m 


(A-6) 


as  follows.  If  we  define  a mode-weighted  average  of  a function  f(z)  by 


(A- 7) 


then  wo  can  compare  the  normalization  for 
6,  Eq.  (3), 

D 


m 


Eq. 


(A-2),  with  that  for 


2 2 

(N  -u  )dz 


1, 


to  derive 


[N2]  - y2([n2]  -j) 

mm  m 
or 

Y2  - (1  - ^/[N2])-1  (A-8) 

m m 

Now,  first  order  perturbation  theory  applied  to  the  self-adjoint  equa- 
tion 


II 


N 


0 


provides  a relationship  between  small  changes  in  the  parameter  k and 
eigenvalue  c (Courant  and  Hilbert  1965,  pp.  353-6): 


[N2]  6c"2  - 6k2  = 0. 
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With  a little  algebra,  this  becomes 
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and  by  (A-6) , 
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Thus,  according  to  (A-8), 
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